Self-energy renormalization around the flux phase in the t — J modeh Possible 

implications in underdoped cuprates 
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The flux phase predicted by the t — J model in the large- A'' limit exhibits features that make 
it a candidate for describing the pseudogap regime of cuprates. However certain properties, as for 
instance the prediction of well defined quasiparticle peaks, speak against this scenario. We have 
addressed these issues by computing self-energy renormahzations in the vicinity to flux phase. The 
calculated spectral functions show features similar to those observed in experiments. At low doping, 
near the flux phase, the spectral functions are anisotropic on the Fermi surface and very incoherent 
near the hot spots. The temperature and doping evolution of self-energy and spectral functions are 
discussed and compared with the experiment. 

PACS numbers: 74.72.-h,71.10.Fd, 71.27.-|-a 
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The pseudogap (PG) phase of high-Tc cuprates 
presents unusual properties in clear contrast to those ex- 
pected for conventional metalg"'^ and, in spite of the many 
scenarios proposed, the problem remains open. A well 
known phenomenological theory is the d-CDW— in which 
a normal state (NS) gap was proposed. Interestingly, d- 
CDW is closely related to the flux phase (FP) which was 
formerly predicted in the context of the t — J modelj^ 
In both scenarios the PG coexists and competes with 
superconductivity, is complex and shows d-wave symme- 
try. FP was also studied by means of different analyt- 
ical and numerical methods.**' In the mean field ap- 
proximation of the t — J model, below a characteristic 
temperature T*, a c?-wave NS gap opens on the Fermi 
surface (FS) near the hot spots (k = (0,7r)). This leads 
to pockets with low spectral weight intensity in the outer 
section.^ Thus, these pockets resemble the arcs observed 
in ARPES experiments,^ Furthermore, recent progress on 
the t ~ J model shows that the competition between FP 
and superconductivity leads to a phase diagram qualita- 
tively similar to that observed in cuprates.^ Raman and 
tunneling spectroscopy results show similarities with the 
experiment 1^ Very recent ARPES Raman specific 
heatfi^ and femtosecond optical pulses^ experiments to- 
gether with some theoretical developments^- support the 
competing two-gap scenario. From the above reasons 
FP may be considered as a good candidate for describ- 
ing the PG region. Nevertheless there are some draw- 
backs to be addressed. While in the FP scenario a true 
phase transition occurs at T* , many experiments show a 
smooth behavior as a function of temperature fi^ In addi- 
tion, well defined quasiparticle (QP) peaks are predicted 
everywhere on the Brilloiun zone (BZ), above and be- 
low T*, while experiments show that the QP exists only 
near nodal directionjil Near hot spots, coherent peaks 
are observed neither below nor above the PG tempera- 
ture in apparent contradiction with the FP picture. On 
the other hand, the PG observed in ARPES, consistently 
with tunneling experiments,*^ seems to be filling up but 
is not closingfi^ leading to a smooth evolution of spectral 



properties. 

In spite of the above objections we will assume the 
point of view that the basic physics of cuprates is con- 
tained in the t — J model (see Ref.pO] for discussions). 
Moreover, to address the above discussion it is neces- 
sary to perform a controllable self-energy and spectral 
function calculation beyond mean field. Recently it was 
proposed a large- approach for the t — J model based 
on the path integral representation for A-operators.-^' At 
mean field level the approach of Ref. 6| agrees with slave 
boson^i^ and the approach of Ref. [Hi- Additionally, it 
can be easily applied beyond mean field allowing the cal- 
culation of self-energy and spectral functions. In this 
paper we discuss spectral functions and self-energy cal- 
culations in the NS in the vicinity to the FP instability. 
We have found that the spectral functions and the cor- 
responding self-energies are very anisotropic on the FS 
leading to well defined QP near nodal direction and very 
incoherent features near hot spots. 

We study the two dimensional t ~ J model with hop- 
ping t and Heisenberg coupling J/t — 0.3 between nearest 
neighbors sites on the square latticci^ In the following 
energy is in units of t. Using the formulation described in 
Ref. Ji] the mean field homogeneous Fermi liquid (HFL) 
becomes unstable when the static {oj = 0) flux suscepti- 
bility Xfiux{ci,u;) = (|)2[(8/J)A2 -n(q,L^)]-i diverges. 
n(q, oj) is an electronic bubble calculated with a form fac- 
tor 7(q, k) = 2A{sin{k^ - - sin{ky - 9y/2)).^>2i.23 
(For definition of A see below). In Fig. la the phase dia- 
gram in the doping-temperature (S—T) plane is presented 
where T* indicates the FP instability onset. At T = a 
phase transition occurs at the critical doping 6^ ^ 0.12 
for the incommensurate critical vector qc = (l,0.86)7r 
near (tt, tt). At finite T the instability is commensu- 
rate (qc — (tt, tt)). Therefore, since the instability takes 
place at, or close to, (tt, tt) the form factor 7(q, k) trans- 
forms into ~ {cos{kx) — cos{ky)) which indicates the 
d-wave character of the FP. Since T* falls abruptly to 
zero at Scr^ we associate the FP with the scenario dis- 
cussed in Fig. la of Ref.jTl]. In Fig. lb we have plotted 
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FIG. 1: (Color online) (a) Large- A'^ t — J model phase diagram 
in the S — T plane, (b) Softening of the d-wave flux mode 
for several 5 values approaching Sc — 0.12. (c) QP weight 
Z for EijA (black circles), T,fiux (red squares) and Et (blue 
triangles), (d) Zfiux versus the FS angle (j) ioi S = 0.13. 



Imxfiux{''\ = (tt, tt), w) at T = which shows that when 
5 ^ 5c & low energy d-wave flux mode becomes soft ac- 
cumulating large spectral weight. This soft mode reaches 
ijj = {) dX 5 = 5c freezing the FP. 

For studying the soft mode influence on one electron 
spectral properties we calculate self-energy renormaliza- 
tions in the NS in the vicinity of T* . The inclusion of 
the superconducting ground state in the calculation does 
not change significantly the FP region (see Refs.jl] and 
[l0|). After collecting all contributions in 0{1/N) the 
self-energy can be written as^^ (in 0{1/N) the inclusion 
of self-energy corrections in the calculation of electronic 
bubbles is not necessary) 

/mI]T(k, w) = ImYjRx{)<i,uj) + ImEfiux{K^) (1) 

where 

ImJ:Rx{k,u;) = -^^{n'^ Im[Dnii{q,uj - ek-q)] 

+ 2n Im[DxR{ci,uj - ek-q)] 
+ Im[Dx\{ci,uj - ek-q)]} 
X [npi-ek-q) + nsiuj - ek-q)] , 



(2) 



and 



X [nF{~ek~q) -f n_B(w - ek-q)] 



In the above expressions, VL — {sk-q + + /i)/2, Ek = 
{t5 + A)(cos(fcj,) -I- cos{ky)) is the mean field electronic 
dispersion, ek ^ Bk - fJ. and, A = J2k cos{kx)nF{,ek). 



ns (np) is the Bose (Fermi) factor, /i is the chemical po- 
tential and Ns is the number of sites. corresponds to 
the charge {5R) sector, nondouble occupancy {5X) sector 
and the mixing of both.^^'^^ Since for J = 0.3 there is no 
important influence of J-contributions in Tirx then, S^ja 
is close to the self-energy for J = 0. Instead, S//„^ is 
dominated by J. The explicit expressions of Dbr, Dxr 
and Dxx are given in Ref . [2l| . 

In Fig.lc we have plotted the QP weight Z = 1/(1 — 
dReY./duj) vs 5 for each self-energy contribution. Zrx 
(black circles) are results at the Fermi crossing momen- 
tum in the (0,0)-(0,7r) direction (antinodal-Zci?). Zrx 
decreases (linearly) with decreasing doping (Zrx — * 
when (5 — > 0). Due to the fact that T,rx is very isotropic 
on the FS, results for Zrx at the Fermi crossing momen- 
tum in the nodal direction (nodal-fc^?) are very similar 
to those for the antinodal-fci^-. In contrast, Zfiux (red 
squares) is strongly anisotropic on the FS. For nodal- A: 
(open red squares) Zfiux 1 which means that ^fiux 
has no important contribution near nodal direction. For 
antinodal-Zcp (closed red squares) Zfiux is rapidly de- 
creasing with doping, i.e., Zfi^x — > when 5 5c where 
the flux instability takes place. In order to show explic- 
itly the anisotropy of Tijiux, in Fig. Id we have plotted 
Zfiux versus the FS angle (j) from ^ = (antinodal-ZcF) 
to (p — 7r/4 {nodal-kp), for 5 = 0.13. The soft d-wave 
mode (Fig. lb), which scatters electrons on the FS with 
momentum transfer q ~ (tt, tt), is responsible for the ob- 
served anisotropy. Using Zrx and Zfiux the total QP 



weight can be written as Zt = ^/{Z 
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Blue triangles are the results for Zt at the antinodal-fci?. 
While for large 5 self-energy renormalizations are domi- 
nated by SjjA, "^fiux dominates near 5c. Results for Zt 
at the nodal- /cf (not shown) follow the same trend as 
Zrx, j-e., the nodal Zt exhibits a linear behavior with 
slope similar to that discussed in Ref. 26]. 

Next we will consider the T,fiux contribution to the 
spectral function AfiuxO^^Lu). Spectral functions for dif- 
ferent 5 values near 5c at the antinodal-fci? are presented 
for T = in Fig. 2a. When 5 ^ 5c, the area under 
the QP peak decreases (consistently with the behavior 
of Zfiux) and spectral weight is transfered to the inco- 
herent structure at about oj ~ 0.05 — O.lt. Fig. 2b shows 
the Lo behavior of the flux scattering rate /mS fiux at the 
antinodal-^F for several 5 values near 5c. With decreas- 
ing doping towards 5c, ImTifiux evolves from a FL-like 
behavior uP') (see solid line for 5 = 0.20) to a non FL- 
like behavior w) (see double-dotted-dashed line for 
5 — 0.13). At the same time, the structures (marked with 
arrows) that appear at low energy are correlated with the 
soft mode discussed in Fig. lb. Im'Efiux at the nodal-kp 
(not shown) is weaker than for the antinodal-Zci? and, in 
addition, it shows a FL-likc behavior for all 5 values. In 
(3) summary, as observed in experiments at low dop- 



ing the scattering rate is strongly anisotropic on the FS 
presenting a FL-like (~ w^) behavior near the nodal-/cF 
and a non FL-like (~ lu) behavior near the antinodal-fci?. 
In addition, this anisotropy disappears with overdoping 
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FIG. 3; (Color online) io = ^ spectral function intensity at 
the antinodal-ZcF in the S — T plane. 



FIG. 2: (Color online) (a) T = spectral functions at the 
antinodal-ZcF for several doping values approaching Sc- A 
broadening rj — 0.01 was used in the calculation. The ver- 
tical line marks the Fermi level, (b) T = flux scatter- 
ing rate for several S values at the antinodal-^F- (c) and 
(d) Spectral functions at the antinodal-^F for T = 0.025 
and T = 0.035 respectively, for 5 — 0.16 (red dashed line), 
(5 = 0.13 (black double-dotted-dashed line) and 5 = 0.10 (vi- 
olet double-dashed-dotted line), (e) The same as in (b) but 
for r = 0.015. 



where the behavior is more FL-Uke. 

In Figs. 2c and 2d spectral functions at the antinodal- 
kp are presented for T — 0.025 and T = 0.035, respec- 
tively. For a given temperature the intensity at the FS in- 
creases with increasing doping showing, consistently with 
experimentSfiii^l a more clearly defined QP peak with 
overdoping. It is also interesting the following behav- 
ior. For instance, for T = 0.035, with decreasing doping 
while the peak at the Fermi level loses intensity a broad 
shoulder appears at positive uj together with an apparent 
gap-like (which is not true gap) structure. As suggested 
by ARPES experiments, ^'"'^^ with increasing temperature 
the gap- like features seem to be filling up but are not clos- 
ing (see results for 6 = 0.10). Notice that for S — 0.10 the 
leading edge below the FS is at w = — 0.05t which, using 
t = OAeV, is about 20 — 30mey thus, of the order of the 
observed PG valucj^^ Therefore, spectral functions arc 
strongly temperature dependent, even at temperatures 
much lower than the hopping i, showing broad features 
at the antinodal-fci?. Spectral functions are nearly doping 
and temperature independent at the nodal-kp and show 
well defined QP peaks as in experiments. ^^'^^ As men- 
tioned above, at mean field level when the temperature 
decreases a gap opens at T* near hot spots. However, 
beyond mean field, spectral functions are very incoher- 
ent above T* showing that features expected to occur for 
T < T* already appear for T > T* due to self-energy 
renormalizations. For instance, the coherence loss when 
going from nodal-kp to antinodal-fci? may suggest the 
appearance of arcs. Therefore, we think that self-energy 



corrections could be masking a likely abrupt change at T* 
given the appearance, in agreement with observations, of 
a smooth crossover. In addition, it may probably be that 
T* determined from the ARPES line-shape is somewhat 
shifted from the temperature at which PG opens. 

The above results are due to the doping and temper- 
ature behavior of the self-energy. In Fig.2e we present 
results for ImYifiux at T = 0.015 for several 5 values at 
the antinodal-fci?. While for 5 = 0.20 Im'Efiux shows, 
as expected for a FL, the maximum at the Fermi level, 
with decreasing doping (see for instance double-dotted- 
dashed line for S = 0.13) the maximum is shifted respect 
to w = 0. As discussed in Refs.,29i,30j. this behavior can 
damage the quasiparticle FL picture. For nodal-fc^ the 
situation is FL-like, e.g., ImYifiux shows, for all doping 
values, the maximum ai lj — Q which increases with T 
as ~ T^. Notice the similarities between our self-energy 
results in Fig.2e and those discussed in Ref . [30l| . 

In Fig. 3 the lo — spectral function intensity at the 
antinodal-fci? is presented in the 5-T plane. The spectral 
intensity distribution suggests a crossover from a strange 
metal at low doping and high temperature to a FL at 
high doping and low temperature. Based on different ar- 
guments, recent works have proposed a similar picture.— 

Before closing let us discuss the influence of on 
previous results. In Fig. 4a self-energy results for S_ra, 
Tjfiux and S^, at the antinodal-Zci;', are presented for 
5 = 0.13 and T = 0.02. ImT^Rx (sohd line) shows a FL- 
like behavior near w = 0. The dotted-dashed line is the 
result for ImTiT (eq.(l)). While at large frequency Et 
is dominated by S/ja, properties near the FS are domi- 
nated by Yifiux- ImYiT at the nodal-fci? is very similar 
to ImYijix (solid line) because S/iux has not important 
influence in the nodal direction. In Fig. 4b, using S^, 
the total spectral functions at the riod&X-kp (solid line) 
and at the antinodal-fci;' (dashed line) are presented for 
low Lo. For the spectral function behavior at large uj see 
Ref. [2l|. Although might modify our previous re- 
sults, it remains the main tendency arising from Yifi^xi 
i.e., the spectral functions are strongly anisotropic on the 
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FIG. 4: (Color online) (a) Scattering rate at the antinodal- 
kp for S = 0.13 and T = 0.02 for E/i„^ (dashed line), E^a 
(solid line) and Et (dotted-dashed line), (b) Total spectral 
function (calculated using Et) for 5 = 0.13 and T = 0.02 at 
the nodal-^F (solid line) and antinodal-^F (dashed line). The 
intensity at the nodal-fc^ was cut for clarity reasons. 



FS, losing intensity and becoming broad when going from 
nodal-fci? to antinodal-fc^. 



Concluding, in the framework of the t — J model we 
have computed self-energy corrections in the vicinity of 
the FP instability. The corresponding spectral functions 
are anisotropic on the FS and very incoherent near hot 
spots thus, they show features similar to those observed 
in experiments. In agreement with recent discussions (see 
Ref. 28] and references their in) our results show two con- 
tributions to the scattering channel. One contribution 
{T,fiux) is strongly anisotropic on the FS, the other one 
(Sii;A) is very isotropic, ^fiux dominates at low dop- 
ing (near 6c) and low energy being relevant for the PG 
properties. In contrast, S^^a dominates at large doping 
and large uj. In Ref. [3^ . by means of transport measure- 
ments, it was found that an isotropic scattering channel 
dominates in overdoped samples which agrees with our 
predictions for S^a- In a- recent work^S it was proposed 
that the high energy features observed in ARPES''^ may 
be due to the large to behavior of E^a- Finally, our re- 
sults indicate that some problems, which can be inter- 
preted against the FP scenario, may be less serious after 
including self-energy renormalizations. 

The author thanks to M. Bejas, A. Foussats, H. Parent, 
A. Muramatsu, and R. Zeyher for valuable discussions. 
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